syms V
Pmax = 530.e03; m = 1435; A = 4; CD = .32; rho=1.2;
Vmax = 200.e03/3600;
double(int(m*V/(2*Pmax/3 - (rho*CD*A*V^3)/2),V,[0,Vmax]))